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Preface 


The  woik  herein  was  performed  in  the  Hydraulics  Laboratory  of  the  U.S. 
Army  Engineer  Waterways  Experiment  Stadrni  (WES)  as  part  of  an 
investigation  into  sediment  disposal  in  die  James  River  for  the  U.S.  Army 
Engineer  District,  Norfolk  (NAO).  This  report  presents  the  results  of  the 
three-dimensional  numerical  modeling  work. 

The  work  was  conducted  from  November  1992  to  April  1993  under  tlw 
direction  of  the  following  personnel:  Messrs.  F.  A.  Herrmann,  Jr.,  Chief  of  the 
Hydraulics  Lrdxuatcny;  R.  A.  Sager,  Assistant  Chief  of  the  Hydraulics 
Laboratory;  W.  H.  McAnally,  Chief  of  the  Estuaries  Division,  Hydraulics 
Laboratory;  D.  R.  Richards,  Chief  of  the  Estuarine  Simulation  Branch, 
Estuaries  Division;  and  Project  Manager  R.  A.  Evans,  Jr.,  Estuarine  Simulation 
Branch. 

Mr.  Evans  prepared  this  report,  and  Messrs.  Richards  and  McAnally,  and 
Dr.  R.  C.  Berger,  E.stuaries  Division,  assisted  in  the  analysis  of  the  results. 

Mr.  Mark  Hudgins,  NAO,  served  as  the  District’s  project  coordinator. 

At  the  time  of  publicaticm  of  this  report.  Director  of  WES  was  Dr.  Robert 
W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contetiK  cf  ihit  report  are  not  lobe  used  for  advertising,  publication, 
or  promoiUmal  purposes.  Citadon  cf  trade  names  does  not  constitute  an 
effieial  endorsement  or  approval  for  the  use  of  such  commercial  products. 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Nm^I  units  of  measurement  used  in  diis  iqx»t  can  be  conveited  to  SI  units 
as  follows: 


1  IlMMpiy 

By 

To  Obtain 

1  *'^**ff  fMt  par  Moond 

0.02631685 

cubic  matars  per  sacond 

1  **** 

0.3046 

malsrs 

1  pounds  (fora»)4aoQnd  par  aqum  foot 

47.38026 

pascals^aoond 

1  pounds  (mass)  par  squaia  foot 

4.882428 

Ufognant  par  aquaia  nnalsr  | 

1  Introduction 


Objective 

Several  dredged  material  disposal  alternatives  ate  being  investigated  by  die 
Norfolk  DiAtict  (CENAO)  for  die  Lower  James  River.  They  include,  but  ate 
not  limited  to:  (1)  ofEshote  disposal,  (2)  creadon  of  wetland  disposal  ishmds, 
and  (3)  diqrosal  into  the  deq>  trough  at  the  entrance  of  die  James  River. 
CENAO  has  a  thorough  knowledge  of  ofGshore  disposal  issues  and  is  cunendy 
conducting  in^iouse  investigatkms  into  the  wetland  island  diqiosal  altemadves. 
Hie  third  option,  deq>  trough  disposal,  is  much  mme  difBcuh  to  evaluate  due 
to  the  duee-dfanenskmal  (3-D)  nature  of  currents  in  die  James  River  entrance. 
The  ohjective  of  dds  study  is  to  determine  if  diqiosal  of  dredged  materials  in 
the  James  River  Deqi  Trough  was  feasible. 


Approach 

Due  to  die  stratified  nature  of  the  lower  James  River  estuary  and  the  need 
to  tqiiesent  bed  vdodty  and  diear  stress,  a  tfaree-dimenskwal  numericd 
modd,  RMA-10,  was  u^  to  predia  die  hydrodynamics.  Aldwugh  die 
original  scope  of  work  required  a  limited  resohidoi  fiume-like  schematic 
model,  a  mote  detailed  fiume-like  sdmmatic  model  whkh  incorporated  mote 
of  the  get^rtqiliical  features  and  a  better  rqnesentation  of  the  tidal  prism  of  die 
James  River  was  used.  WhUe  this  modd  will  give  more  accurate  results  dum 
the  limited  resolution  fiume-like  schematic  modd,  it  remains  a  fiume-like 
schematic  modd  and  die  results  should  be  judged  as  such. 

Severd  numeticd  modd  meshes  were  constructed  to  test  different 
conditions.  A  Base  geometry  was  made  whidi  cmisisted  of  prototype  dzes 
and  dqitiis,  with  mote  detafl  in  die  area  of  the  Deep  Trough.  Hgure  1  shows 
die  modd  boundary  siqietimposed  on  die  James  River  region.  Note  die 
location  of  the  De^  T^gh  at  die  entrance.  Natural  botttan  devations  in  die 
Deep  Tiou^  exceed  -85  feet  MLLW.  Three  Plan  geometries  were  constructed 
which  consisted  of  dunging  die  bonmn  efevaticms  in  die  Deep  Trou^  to 
simolate  different  levds  of  dredge  materid  fill  The  Flan  geometry  depdu  ate 
shown  in  Table  1. 
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2  Hydrodynamic  Model 
Validation 


Udal  stations,  major  tributaries  and  major  freshwater  inflows  are  shown  in 
Figure  1 .  Note  that  the  nxxlel  has  two  tidal  boundaries,  one  in  the  Chesapeake 
Bay,  the  other  in  the  Atlantic  Ocean.  The  elements  of  the  model  were  divided 
into  IS  materia]  types  based  on  depth,  locatirm,  and  whether  the  elements  were 
two-  or  three-dimensional.  Table  2  lists  the  hydraulic  material  types  and 
associated  hydraulic  parameters.  Note  that  material  types  6,  1 1,  and  12  are 
listed  as  boundary  elements.  Hydraulic  constants  for  diese  types  were  adjusted 
to  control  the  instability  near  the  boundaries.  Elements  with  no  3-D  layers  are 

2- D.  The  number  of  layers  corresponds  to  the  number  of  elements  in  the 
vertical.  Since  the  boundary  conditions  for  these  boundaries  were  not  known 
and  were  difficult  to  define,  large  eddy  viscosity  values  had  to  be  used.  These 
values  are  numerically  rather  than  physically  based  but  will  not  compromise 
results  near  the  deep  trough  due  to  the  large  distance  from  the  boundaries.  It 
should  be  noted  that  the  Marming's  n  values  are  generally  higher  than  one 
would  expect  from  previous  2-D  modeling  experience.  This  is  a  result  of  the 

3- D  hydr^ynamic  code  using  the  near  bed  velocity  to  calculate  shear  stress, 
not  a  vertically  averaged  velocity. 

Figure  2  shows  the  area  of  interest  (the  Deep  Trough)  with  locations  of 
velocity  stations.  The  stations  JG-xx-yy  correspond  approximately  to  velocity 
stati<Mis  on  the  Cliesapeake  Bay  Hiysical  Model,  The  stations  Tl,  T2,  T3,  TS, 
and  TN  are  in  locations  in  or  near  the  Deep  Trough. 

The  3-D  numerical  model,  RMA-10,  was  run  for  a  total  of  24.3  tidal  cycles 
vtitii  a  repetitive  Mj  tide  (period  =  12.42  hours)  for  a  total  of  303  hours  (12.74 
days).  Using  the  Base  geometry,  a  single  simulation  of  253  hours  was  made 
as  a  spin-up  to  create  hotstart  conditions  for  the  tiKxlel.  This  hotstart  was  used 
with  the  Base  and  Plan  geometries  to  generate  SO  hours  of  simulation.  Only 
the  last  25  hours  of  these  results  were  used  in  the  validation  and  plan 
comparison  in  order  to  minimize  the  effects  of  the  hotstart  conditions.  The 
tidal  amplitude  for  the  last  50  hours  of  the  simulation  was  4.2  feet.  The 
salinity  coiKentrations  at  the  Chesapeake  Bay  and  Atlantic  Ocean  boundaries 
were  22  and  25  ppt,  respectively.  Freshwater  inflow  at  the  upper  end  of  the 
James  River  was  set  at  9,181  cfs.  These  boundary  conditions  are  greatly 
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Table  2 

Hydraulic  Oonatanta 
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1  1 

simplified  representations  of  real  condititHis  that  were  necessary  to  conduct  this 
limited  scope  study. 

Since  this  is  primarily  a  flume-like  schematic  model,  the  numerical  results 
in  the  vicinity  of  the  Deep  Trough  were  the  only  results  considered  in  the 
validation  of  the  model.  Figure  3  shows  the  numerical  model  predicted  tide 
for  Old  Pt.  Comfort  and  Newport  News,  with  the  physical  model  tide  ranges. 

A  comparison  of  tide  ranges  between  the  numeiicd  and  physical  models  shows 
a  difference  of  approximately  0.12  feet  for  Old  Pt.  Comfort  and  0.30  feet  for 
Newport  News. 

Figures  4  through  6  show  the  numerical  nKxJel  surface  and  bottom 
velocities  for  the  three  physical  model  stations  nearest  to  the  Deep  Trough. 

The  dashed  lines  indicate  the  maximum  flood  and  ebb  velocities  measured  in 
the  physical  model  for  the  same  locations.  The  numerical  model  results  are  in 
fair  agreement  with  the  physical  model  results,  with  the  best  agreement  being 
in  the  Deep  Trough  at  station  JG-01-03.  Figure  7  shows  the  tide  ranges 
predicted  by  the  numerical  model  and  compared  to  the  physical  model  for 
three  stations  in  the  upper  pait  of  the  model.  Overall,  velocities  in  the 
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numericaJ  model  are  too  small  compared  to  the  physical  model.  There  are  two 
reasons  for  this.  First,  due  to  the  schematization  of  the  James,  large  amounts 
of  intertidal  areas  were  not  included,  resulting  in  a  smaller  tidal  prism  and, 
therefore,  smaller  velocities.  Second,  vertical  mixing  appeared  to  be  greater  in 
the  numerical  model,  reducing  surface  currents  which  would  have  been 
affected  by  a  larger  salinity  gradient.  These  are  somewhat  related,  but  are  not 
considered  fatal  flaws  due  to  the  limited  expectations  for  a  tiKxleling  effort  of 
this  level. 

Figures  8  through  1 1  present  the  surface  and  bottom  velocities  for  the  Base 
and  Plan  geometries  at  three  stations  in  the  vicinity  of  Old  Pt.  Comfort. 

Figures  12  through  IS  present  the  surface  and  bottom  velocities  for  the  Base 
and  Plan  geometries  at  three  stations  in  the  vicinity  of  Newport  News. 

Figures  1 6  UiTOugh  23  show  the  surface  and  bottom  velocities  at  locations 
in  the  trough  (T1 ,  T2.  and  T3)  and  both  north  and  south  of  the  trough  (TN  and 
TS,  respectively).  These  also  show  very  small  changes  in  the  velocities  due  to 
the  various  plans. 
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3  Shear  Stress  Analysis 


Shear  stress  along  the  bottom  was  computed  using  die  following  ftnmula: 

T  -  1/2 


with 

X  s  shear  stress,  IhAq  ft 
p  B  water  density,  1.935  slugs 
^  =  current  friction  factor,  0.007 
u  s  water  velocity,  ft/sec 

Based  on  information  provided  by  die  Norfolk  District,  typical  sediment 
size  to  be  disposed  in  the  trou^  will  be  qiptDximatdy  0.2  mm  (fine  sand). 
The  critical  shear  stress  (‘O  for  tpiartz  sifo^  of  diis  size  is  approximately 
0.00S  Ib^  ft  This  study  examined  scour  patterns  for  two  values  of  cridr^ 
shear  stress.  Areas  where  die  critical  shear  stress  equalled  or  exceeded 
0.008  Ih^  ft  at  any  time  during  die  last  two  tidal  cycles  woe  examined.  This 
value  is  very  conservative  and  erosion  is  most  jHobi^e  for  sudi  a  hi^  value. 
However,  tl^  value  may  only  occur  at  maximum  ebb  and  flood  and  may  not 
be  rqnesentative  of  a  significant  portion  of  die  tidal  cyde.  In  fact  tills  value 
was  not  exoeetted  for  a  significard  portion  of  die  tidal  cycle.  The  second  value 
of  critical  shear  stress,  0.006  Ifa/sq  ft  was  used  to  develop  scour  patterns.  But 
diese  patterns  repesent  areas  iiriiere  die  shear  stress  exceeded  0.006  Ib/sq  ft  at 
least  33  percent  of  die  time.  While  0.006  Ib/sq  ft  is  a  less  conservative 
estimate  of  critical  shear  stress  dian  0.008  It/sq  ft,  it  still  exceeds  the  value  of 
0.005  OVsq  ft  given  for  quartz  spheres. 

Hgures  24  duough  27  show  regions  widi  shear  stress  greater  than  or  equal 
to  0.008  Ib^  ft  during  some  poim  in  die  last  two  tidal  cydes  of  die 
simulatioiis  for  die  Base  and  Flan  geometries.  Note  that  all  show  regions  of 
potential  scour  coinddem  with  the  Deq>  Trou^  Figure  28  shows  Plan  1 
results  siqimimposed  on  the  Base  results.  Since  Plan  1  is  the  most  drastic 
change  in  geometry,  it  should  show  the  most  deviation  fiom  the  Base.  This 
figure  shows  very  little  change. 

Rgures  24  duougfa  28  show  only  dud.  at  some  times  in  the  tidal  period,  the 
shear  stress  will  be  large  enough  for  scour  to  occur.  However,  they  give  no 
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indication  of  the  amount  of  time  that  critical  shear  stress  will  be  equalled  or 
exceeded.  Figures  29  duoush  33  show  the  regions  in  which  shear  stress 
exceetb  0.006  Ib/sq  ft  at  least  33  percent  of  dw  dme.  It  is  believed  diat  these 
give  a  better  estimate  of  scour  potential.  However,  as  shown  in  Hgures  24 
through  28,  dieie  appears  to  be  ix)  significant  change  in  the  scour  patterns. 
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4  Conclusion 


The  3-diineiidonal  hydrodynamic  model  RMA-10  and  a  schematized  model 
mesh,  widi  average  boundary  condititms,  was  used  to  simulate  water  vdocides 
and  elevadtms  of  the  lower  James  Riva.  It  is  concluded  that  dredge  material 
diqxMed  in  die  Deqi  Trou^  is  not  likely  to  remain  without  protective 
aimming.  This  detorminadon  was  made  by  calculadng  the  bottom  shear  stress 
using  die  near-bottom  water  velocities  and  comparing  widi  a  conservative 
critical  shear  stress  value.  For  eadi  plan  tested,  water  velocities  were 
sufficient  to  cause  scouting  during  parts  of  the  tidal  cyde.  However,  if  a 
mildly  dispersive  site  was  desired,  these  results  indicate  ttiat  it  could  be 
feasible  to  use  the  Deep  Trough  as  a  disposal  area  witixnit  protective  armoring, 
bfore  accurate  results  will  require  a  more  detailed  3-D  modd  mesh  and, 
perhqis,  a  3-D  sedimentatitm  model 
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Figure  4.  Numerical  model  surface  and  bottom  velocities  versus  the  physical  model  veiodty 
ranges  (M,  tide),  Station  JQ^}1-01 
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Figure  5.  Numericai  model  surface  and  bottom  velocities  versus  the  physical  model  velocify 
ranges  (M,  tide).  Station  J&01-02 


Figure?.  Numerical  vereue  physical  model  tidal  ranges.  Upper  James  River 
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Figure  8.  Pt  Comfort  surface  and  bottom  velocities,  Base 
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Figure  10.  R.  Comfort  surface  and  bottom  veioctties,  Ran  2 
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Figure  12.  Newport  News  surface  atxi  bottom  velocities.  Base 
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Figure  13.  Newport  News  surface  and  bottom  velocities.  Plan  1 


Figure  14.  Newport  News  surface  arxf  bottom  velocities,  Plan  2 
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Figure  16.  Base  condition  surface  and  bottom  veiocities,  Stations  T1-T3 
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Figure  17.  Base  condition  surface  and  bottom  veiodties,  Stations  TN  and  TS 
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Figure  18.  Plan  1  surface  and  bottom  veiocWes,  Stations  T1-T3 
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Figure  19.  Plan  l  surtaoe  and  bonom  velodlies.  Stations  TN  and  TS 
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Figure  20.  Plan  2  surface  and  bottom  velocities,  Stations  T1-T3 


Figure  21.  Plan  2  surface  and  txMom  veiocfties,  Stadone  TN  and  TS 


Figure  24.  Base  geometry  -  areas  where  x  2  0.006  lysq  ft  during  the  tidal  cycle 


Plan  1  geometry  -  areas  where  t  i  0.006  fc/sq  ft  during  the  tidal  cycle 


Figure  26.  Plan  2  geometry  -  areas  where  x  ^  0.008  b/sq  ft  during  the  tidal  cyde 


Figure  27.  Plan  3  geometry  •  areas  where  x  2;  0.006  IVsq  ft  during  the  tidal  cycle 


Figure  28.  Base  versus  Plan  1  geometry  -  areas  where  t  ^  0.008  fc/aq  ft  during  the  tidal  cyde 


1.  Base  geometry  region  which  soours  at  least  33  percent  of  the  time 


Figure  30.  Plan  1  geometry  region  which  scoure  at  least  33  percent  of  the  tfene 


Figure  31.  Plan  2  geometry  region  which  sooure  at  least 


Figure  32.  Plan  3  geometry  region  which  scours  at  least  33  percenl  of  the  time 
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